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ABSTRACT 

The objective of this study was to examine the relation- 
ship between forward ship motion and the effects of aero- 
dynamic and hydrodynamic disturbances on motions in the 
horizontal plane. It includes course control and stability 
analysis for the unsteered and steered cases (manual and 
automatic) using several sets of operating conditions. 

Hydrodynamic and aerodynamic effects are considered to 
be functions of hull motion, rudder (steered cases) , pro- 
peller and wind effects ratio. A dimensionless mathematical 
model has been developed and solved with respect to ship 
axes . 
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PHYSICAL QUANTITIES 



Coefficients of perturbation equation in yaw 

Coefficients of perturbation equation in sway 

Coefficients of perturbation equation in surge 

Force on the body 

J Moment of inertia referred to z-axis 
z 

i Ship length 

M Moment about origin 

m Mass of ship 

N Hydrodynamic and aerodynamic yaw moment 

Aerodynamic yaw moment 

N^, Derivative of hydrodynamic yaw moment with respect to 
yaw acceleration 

N^ Derivative of hydrodynamic yaw moment with respect to 
sideslip velocity 

Ng Derivative of hydrodynamic yaw moment with respect to 
rudder angle 

R Vector distance of a point from the origin, and the 

components along the X, Y, and Z axis R = ix + jy + kz 
r Yaw rate 

U Velocity of the origin U = lu + 3^v + Jew 
u^ Component of equilibrium ship speed in x-axis direction 

u Component of ship velocity relative to air in x-axis 

a. 

direction 

U_ Ship velocity relative to air 

cL 

Wind velocity relative to earth 
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V 



Component of ship speed in y-axis direction 

Y Component of ship velocity relative to air in y-axis 

Si 

direction 

v^ Component of equilibrium ship speed in y-axis direction 

X Hydrodynamic and aerodynamic force component in x-axis 
direction 

X Aerodynamic force component in x-axis direction 

Si 

Xp Hydrodynamic force component in x-axis direction due to 
propeller 

X^ Derivative of hydrodynamic force component in x-axis 
direction, with respect to surge acceleration 

X^^ Second derivative of hydrodynamic force component in 

x-axis direction, with respect to sideslip velocity and 
yaw angular velocity 

X^ Drag coefficient 

Xgg Second derivative of hydrodynamic force component in 
x-axis direction, with respect to rudder angle 

Y Hydrodynamic and aerodynamic force component in y-axis 
direction 

Y Aerodynamic force component in y-axis direction 

Y^ Derivative of hydrodynamic force component in y-axis, 
with respect to yaw rate 

Y^ Derivative of hydrodynamic force component in y-axis 
direction, with respect to sideslip acceleration 

Y^ Derivative of hydrodynamic force component in y-axis 
direction, with respect to rudder angle 
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S Rudder angle 

Equilibrium rudder angle 
e Density of water 

e Density of air 

ijj Heading angle of ship 

Direction of wind velocity relative to ship 
Direction of wind velocity relative to earth 
n Angular velocity of the body about the origin 

n.= ip + jq + kr 
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DIMENSIONLESS FORMS 



Quantities 



Typical Symbol 



Typical 

Dimensionless Form 



Length 


X 

r 


X' 

r 




X /n 
r' 


Force 


Y 


Y* 


- 




Moment 


N 


N’ 


= 




Mass 


m 


m' 




m/^.3 


Angular velocity 


r 


r’ 
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rZ/U 


Static force rate 


Y 

V 


Y' 

V 
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Y /®)l 

V 2 


Static moment rate 


N 

V 


N’ 

V 


= 




Rudder force rate 


Y 

s 


Y’ 

s 
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Damping force rate 


Y 

r 


Y* 

r 
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Damping moment rate 


N 

r 


Nr 


= 


Vl‘ 


Inertial coefficient 


Y* 

V 


I 

Y* 

V 
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Moment of inertia 


I 

z 


I* 

z 


= 




Acceleration 


« 

u 


u' 


= 


UJl/U^ 


Velocity 


u 


u* 


= 


U/U 


Angular acceleration 


• 

r 


• 

r' 


= 


rt^/u 


Density of air 


‘'a 




= 




Inertial coefficient 


N‘ 

r 


1 

N* 

r 


= 
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I. 



INTRODUCTION 



In the last twenty years many mathematical derivations 
have been performed in the area of motion of vehicles in a 
fluid medium. Some of these mathematical models have been 
derived using the six degree of freedom equations of motion. 

In this study of six degree freedom equations are re- 
duced to three degrees of freedom with several assumptions. 

Low ship speed is assumed and coupling effects of roll, 

i 

pitch and heave motion on to the horizontal plane motion 
are assumed negligible. It is also assumed that atmospheric 
motion is uniform. 

Aspects of the problem are simultaneous consideration of 
hydrodynamic and aerodynamic effects as well as nonlinearity 
in the external disturbances (especially in aerodynamic terms) . 
Equilibrium conditions and perturbation equations are derived 
from the mathematical model. Nonlinear parts of the mathe- 
matical model are linearized with respect to equilibrium 
conditions . 

The effects of wind disturbances on ship motions are 
studied. A digital computer (IBM model 360/67) was used to 
solve the transient response in yaw and sway for three dif- 
ferent wind directions (head, beam, stern) and constant wind 
velocities with constant ship speed. Also the eigen value 
theorem was used to study ship stability. All these com- 
putations were made for the unsteered (zero degree rudder 
deflection) and steered case (with manual as well as auto- 
matic control) . 
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When wind velocities are several times the ship speed, 
even moderate winds create difficulty in controlling the 
ship course and stability (at low ship speeds) . These ef- 
fects may or may not be controllable, since the hull hydro- 
dynamic reactions are reduced and ship motion is greatly 
effected by wind force, propeller thrust and coupled rudder 
force. This situation can cause unusually large hull drift 
angles extending to between plus or minus 180 degrees. The 
usual hull hydrodynamic force representation is not applica- 
ble beyond 15 to 20 degrees. 
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li. THE MATHEMATICAL MODEL 

A. COORDINATE SYSTEMS AND SHIP MOTION 

The coordinate system is defined with respect to the 
ship itself. It's coordinate origin is chosen to be the 
center of gravity of the ship as shown in Figure 1. 

Positive X axis (lengitudinal axis) is directed to the 
forward direction along the center line of the ship. 

Positive Y axis (transverse axis) is directed to star- 
board, perpendicular to the plane of symmetry. 

Positive Z axis is directed downward, and lies in the 
vertical symmetry plane, perpendicular to the X-Y plane. 

Yaw is a retatienal motion of the ship about the Z axis. 
(Assiime that the ship is moving in the positive X direction 
and that a moment is applied about the Z axis tending to 
rotate or turn the ship in clockwise direction as indicated 
in Figure 2 which shows the positive yaw) . 

1 . Equations of Motion (on the horizontal plane) 

The dynamic response in six degrees freedom for the 
motion of a rigid body [7] is given by equation (1) . 

X = W[U-RV+QW-X^ (R^+Q^)+Y^ (PQ-R)+Z^ (PR+Q) ] 

b bo 

Y = W[V-PW+RU+X^ (R+PW)-Y^ (P^+R^)+Z, (RQ-P) ] 

b b b 

Z = M[W-QU+PV+Xg (PR-Q)+Yg (P+QR)-Zg (Q^+P^) ] 

L = PIj^+(l2“Iy)QR+PI[Yg (W-QV+PV)-Zg (V-PW+RU) ] 

M = QIy+(Ij^-l2)PR+W[Zg (U-RV-QW)-Xg (W-QU+PV) ] 

N = Rl2+(lY“Ij^)PQ+P5[Xg (V-PWtRU)-Yg (U-RV+QW) ] (1) 
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These equations represent the reaction of the rigid 
body. They do not include external forces such as rudder 
force, wind force, aerodynamic and hydrodynamic force of 
water etc. Applied external hydrodynamic and aerodynamic 
forces are considered to be functions of hull motion, rudder 
angle, propeller and wind. 

Assuming that roll, pitch, and heave motions may be 
neglected on the horizontal plane and that x^, y^, z^ are 

D D D 

zero for. the center of gravity at the origin of the system 
and noting that on the horizontal plane U = lu + 3^v (linear 
velocity) , then equation (1) becomes 

X = w[u - Rvl (Surge) 

Y = M[V + RUl (Sway) 

Z = 0 

L = 0 

M = 0 

N = RI^ (Yaw) (2) 

These equations are nonlinear. Analytic solution and 
computer simulation are very difficult. But all nonlinear 
terms might be linearized by using the Taylor series expansion. 

Briefly, small perturbations are considered for: 

X = m (u - rv) 

Y = m(v + ru) 

N = rl^ (3) 



16 












1 ^ 



# 











I 



« 




# 



The external force and moments to be included are; [4] 



X = X(hull) + X (rudder and prop.) + X(wind) 

Y = Y(hull) + Y (rudder and prop.) + Y (wind) 

N = N(hull) + N (rudder and prop.) + N(wind) (4) 

So N, X, Y represent total hydrodynamic and aerodynamic forces 
and moments which are generated by the ship motion, rudder, 
propeller and wind. The following assumptions permit elimina- 
tion of relatively unimportant terms in the hydrodynamic and 
aerodynamic expressions: [1] 

1. Ship asymmetry with respect to the centerline plane 
is minor (this permits elimination of odd ordered terms for 
the X force expansion and even ordered terms from the Y and 

N expansion) . 

2. All terms higher than third order are unimportant. 

3. Rudder force and moment derivatives of higher than 
first order, and effects of rudder angular rate are negligible. 

4 . The lengitudinal component of ship speed is much 
greater than the lateral component (u>>v) . 

5. Based on experimental and analytical results, non- 
linear hydrodynamic terms in the yaw and sway equation for the 
stability analysis are considered minor, except for certain 
nonlinear hydrodynamic terms in the surge equations such as 

X ' . 

vr 

6. Atmospheric motion is uniform. 

Thus equation (4) may be written as: [1] 

X = C,u+C.u^+C,vr+C.u^5^+x,v U +x 
12 3 4 laap 
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y = B, v+B^uv+B-,ur+b..U'^ 6+Y, v_U 

2 






N = Aj^r+ A 2 UV+ A3ur+ A^v'^ 6+N3 U^v^+N2V^U^ (5) 

All hydrodynamic and aerodynamic constants are defined in 
the appendis. 

The hydrodynamic effects, forces are defined as a 

function of the acceleration and velocities of the hull. 

The rudder deflection will also produce an additional 

2 

resistance, an axial force proportional to S . This force is 
expected to be small. The lateral force and moment are pro- 
protional to S. [5] 

Forces and moments due to wind acting on the above 
water protions of the hull act on the ship in the horizontal 
plane. It is necessary to determine the relative wind speed 
and direction with respect to the ship in it's own reference 
frame . 

Magnitude and direction of wind is shown on Figure 1. 
The relative wind speed components along the ship axes are 



u = u+U, cos (t{) +!{;) 
a a. a 



V = v-U sin (4<^+ii)) 

ct d d 



Then 



2 . __2,% . _^^.,-J-'a 



U = (u + V ) ^ , it = tan 
a a a' ' ^a 



in X direction 
in Y direction 

-I'^c 



(6) 



u. 



U and ijj are wind velocity and direction along the earth 
a a 



axis. 
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Figure 1. Relation of Wind Velocity and Direction to the 
Ship Axes. 
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From equations (3), (4), and (5), accelerations in 

yaw, sway, and surge are determined: 

So 

X = u(m-xu) - mrv 
y = V (m-yv) + mru 

N = r(I^ - N^) (7) 

The relationship between velocities of the ship with 
respect to the earth axis and velocities of the ship with 
respect to the body axis is: [3] 

Y^(t) = u sin (i|;) + v cos (\p) 

X^(t) = u cos (i|j) - V sin (4)) 

also 

4-(t) = r (8) 

Where is the heading angle of the ship relative to the 
original path. 



origincj-V 





Figure 2. Positive Direction of Yaw. 
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( 9 ) 



= /(u.sin (i;;) + v.cos (ifi)) dt + Y^ (o) 

= /(u.cos (4)) - V sin (if))) dt + X^(o) 
when if) is small; 

sin (4))=if) / cos (if))=l U=Ue (equilibriiom case) (10) 
So 

Y^(t) = Ue • if; + V 

(t) = Ue - v.if; (11) 

Combining equation (5) and (7) and putting in nondimensional 
form. [2] 

(m'-xu' ) [ (x'+x '6^)u^ + il(x '+m' ) vr+p ’x'U^+x'U^] 

^ (m' -y • ' ) v=^l^ [y ' uv+jj, (y ' -m' ) ur+y ' 6u^+e ' y 'U^ ] 

^ VzV IT 333 

(I'-N* ' )r=^«.3 [N'uv+AN'ur+N' 5u2+e'N'u2 (12) 

z zr z V r 333 

where 



r^£3u2N'=N =5-S,3 u2n 

2 a a a 2 



I 

3 



N' N' 

N' = P'U'2 (^^^sin(2if; )+_^2gin(^ ) ) 

3 3 3 Z 3 Z 3 



^^2u|y;=y^=^i2u2y; 



y' = plU'^Y' sin(if; ) 
■'a ^a a ai a 



^)12u2x'=x =£ji2u2x' 

2 a a a 2 a 



X' = p'U'2x' cos (if; ) 
a *^a a a^ ^a 



( 13 ) 
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u=U since u>>v) 



p; = P^/P , 

Expanding in Taylor's series around the equilibrium conditions 



u = ue 



V = ve 



= 6e 



r = x = x = u = v= =0 (14) 

Substitute equation (14) into equation (12) 

N'ueve + N' 6e ue^ + N'ue^ = 0 
V (S a 

y'ueve + yl6e ue^ + y' ue^ = 0 
■*v ' ■'6 ■''a 

x'ue^ + x.-'Se^ ue^ + x' ue^ + x'ue^ = 0 (15) 

o 66 a p 

Sideslip velocity and rudder angle in the equilibrium condi- 
tion can be expressed as: [1] 



Ve' 





Y 

V 



^6 



I 



6e 



N 'Y ' 
a V 



N 'Y.' 
V 6 



Y 'N ' 

a V 



Y 'N^ ' 

V 6 



(16) 



Perturbation terms (small deviation from the initial straight 
course) may be denoted by bars. So actual motion is 

ue = ue + U V = ve+V S = Se + S 

= u = u v = v (17) 

Perturbations in aerodynamic quantities due to , v' and u' 
are ; 
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N' 

a 



Y' 

a 



3N' 3N' 

a 7 

4 ) + 

3 i|; 3 v ' 



^ V- + 



3Y' 



3.Y' 

a 7 , 

^ + 

3 v * 



a ~ I . 

— V* + 



3N* 
3 u' 
3 Y' 
3u ’ 






a u« 






sx' 

i|) + 

34 » 



3X' 



— v' + 



3X' 



a - , 
— u' 



3v' 3u' 



( 18 ) 



Where 



3N' 

3. 

3 i|) 



= p'[(N' cos 2^b + N' cos 4) ) (U'2 - u') 

33]^ 3 32 Si 3 3 



N' 

-2U' sin 4» ( — ^ sin 24) + N' sin 4)_)] 



9N' 
a 

3v' 



N' 



= -p' U ' [ (N ' cos 24) + N' cos 4) ) + 2 (— ^ sin 2\p 



3 3 3 2 



az 



+ N' sin 4)^) sin 4) ] 
32 ^ 3 - 



3N' 

— — = p' U' [ (N' COS 2ip + N' cos 4) ) (-sin 4) ) 

3 a aj a az a a 

N' 

+ 2 (— — ^ sin 2i> + N' sin 4» ) cos jp ] 

^ 3 32 3. 3 



3Y' 

— ~ ~ ” 2 U' sin^ 4) ] 

3 3^ 3 3 3 3 3 



3Y' 
a 

3v' 



■‘"i “A ♦a> 



3Y' 
a 

3u' 



‘^a “i *a *a 
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ax' 

a 

ax 



- p ' X' sin (U' 2 + U' ) 

3. d. 1 3. 3i 3 



av' 



- -o' X' u' cos 4/ sin 4/ 
*^a a^ a a ’^a 



ax' 

a 

au' 



= p' X' U’ (1 + COs2 4; ) 

a. 3. a. 3 



( 19 ) 



The coefficients in equation (19) where taken from reference 
[2]. Aerodynamic coefficients are shown in Figure 3. 

Substitute equation (17) and (18) into equation (12) 
and obtain the perturbation equation in the following form; 



• 2x' 2x' 2x ' 

(m'-x* ' )u'=(6 '+m' )v^'r'+ — — ^v' + (2X ' +2X» — - 

u vr e 2x 2v' ° ® 2u' 



) 



G'+2Xg^'6^6 



2Y' 2Y' 2Y' 

(m ' ' ) V ' = (y^ ' -m ' ) ? ' +— (y^^ ' +—2.) v ' + (y ’ v ' +2y ’ 6 e+— 2.) 

^ 2x ^ 2v' ® 2u' 



u'+y^6 



• 2N' 2N' 2N' 

(I'-N* ' )r'=N'r’+ — ^ij;+(N'+ — ^)v'+ (N' v'+2N;6 + — ^) 
^ ^ 24 ^ 2v' ® ® 2u' 



u'+NlS 

0 



( 20 ) 



The mathematical model which is shown in the above equations 
is developed only for the forward ship motion. The stern 
ship motion equations are v/ritten by putting negative signs 
in front of some of the coefficients. 
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MVOHOOYNAMIC CO£FFIClENT3 




M-60R2 M-60 M-60 R I 

RUDDER FORCE COEFFICIENT 



Fxgure 3. Hydrodynamic Coefficients Vs Rudder Force Coefficient. 
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Example: (derivative of hydrodynamic yaw moment with 

respect to sideslip velocity) takes a negative value for the 
forward motion, positive value for the stern motion. 

(derivative of hydrodynamic force component in 
Y-axis, with respect to yaw rate) takes negative value for 
the forward motion, positive for the stern motion but is not 
a sensitive parameter since m' is the predominant term in 
(Yj.' - m') . 



26 









■m $ 



f 







III. CONTROL STUDIES 



The linearized equation of motion was derived and found 
in equation (20) . In this part two different problems were 
simulated. One of them was manually steered ship motion in 
wind disturbance. The other was automatically steered ship 
motion under the same effect. So computer program (II) and 
(III) were developed for these purposes respectively. The 
computer model is shown in Figure 4, for these programs. 
Program (II) and (III) are similar to each other. There is 
only one difference in section 4 of both programs. Program 
(II) was integrated using the second order equation, however 
program (III) was integrated with the third order equation. 
This difference was caused by the feedback loop (page 38 ). 
Also the Laplace Transform Technique could have been used to 
obtain the system description for simplicity and system 
characteristic equation for the parameter plane program. The 
characteristic equation was calculated in the computer pro- 
gram. Also three more additional and different mathematical 
sections were calqualted in the same program. One of these 
sections included the parameter calculation from equation 
(6) through equation (20) . Secondly the transfer functions 
of yaw and sway were calculated from the three degree of 
freedom equation. The last section was the polynomial cal- 
culation for the eigen values. 

The last program, (parameter plane (IV) , was used to get 
optimal values for yaw gain constant and yaw rate gain con- 
stant using the constant zeta (function omega) and constant 
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omega (function zeta) curves for automatically steered ship 
motion. According to reference [1] and [6] , these optimal 
values were chosen from Figures 40 through 44 assuming 
zeta =• 6, omega = 11.5 for varying wind direction and 
velocities. 

The transform of the linearized equations of motion (20) 
from the T domain to the S domain, for all initial conditions 
equal to zero are: 

(ABC.<S2-AAD.iS-BFA)ij;-(BGD.<S')- -(BGC.5)- = AAE.6 

6 6 

(-BGE.S-BFD)iJJ- (BGG.S)- + (ABD .^2 _bgf . 5 ) - = AAC . 6 (21) 

6 6 

(-CAA.-S-BFG)iji+(ABE.S2-CAB.S)- -(BGA..S)- = CAC . 6 

6 6 

In equation (21) different letters are used for simplicity. 
For example: ABC = 1' - N* ' , AAD = N' , BFA = , 

ABE = m' ” , etc. 

let 




B 

Applying Cramer's rule, the transfer function for — and— 

S S 

are : 





AAE 

AAC 

CAC 


-BGD . S 
-BGG.S 

ABE.S2-CAB.S 


AAE 

AAC 

CAC 


s 


( ABC . S 2 -AAD . S -BFA ) 


-BGD.S 


-BGC . S 




(-BGE.S-BFD) 


-BGG.S 


ABD.S2-BGF.S 




(-CAA.S-BFG) 


ABE.S2-CAB.S 


-BGA.s 
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The open loop transfer function is; 



ijj WA . S2 +WC . S+WF 



6 ssa.s‘»+ssb.s3+sse.s2+sbs.s+ses 
W here 

WA = - ABC.ABD.ABE etc. 

All terms were calculated using program 1. 

In this case, the results are a second order niimerator 

and fourth order denominator. 

B(s) 






Transfer function 



S(s) 



, was also obtained using Cramer's 



rule. 



B 

I 



ABC.S^-AAD.S-BFA 


-BGD.S 


AAE 


-BGE.S-BFD 


-BGG.S 


AAC 


-CAA.S-BFG 


ABE. S^ -CAB. S 


CAC 


abc.s2-aad.s-bfa 


-BGD.S 


-BGC.S 


-BGE.S-BFD 


-BGG.S 


ABD.S^-BGF.S 


-CAA.S-BFG 


ABE.S^-CAB.S 


-BGA.S 



(23) 



The transfer function between B and S is: 



B WWA.S^+WWC.S^+WWF.S+WVB 

- (24) 

6 SSA. S^+SSB . S‘^+SSE . S 3+SBS . S^+SES . S 

where 

WWA = -ABC.AAC.ABE etc. 

Again all terms were calculated using program 1. 
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For the lateral displacement, sway is defined by 




The transfer function is 

Y WWA.S +(\'7WC+WA)S + (V7WF+WC) S+WVB+WF 



S SSA.S +SSB.S +(SSE+SSA)S +(SBS+SSB)S +(SES+SSE)S +SBS.S+SES.S 

(25) 



A. MANUALLY STEERED SHIP MOTION 

As mentioned before, the goal of this section is to get 
the ship trajectories in constant wind effect with low ship 
speed. For this purpose computer program (II) was developed 
and applied to equation (20) . Two situations were considered; 
Situation 1; zero degree rudder deflection (unsteered). 
Situation 2: positive 15 degree rudder deflection. 

All these situation were studied assuming that the ship moved 
in the forward direction at a 1 mile/hr speed in. a calm sea. 
Each of the wind velocities was applied suddenly. 

Computer program II responses are shown in Figures 5 
through 30 for yaw and sway. These results were also tabu- 
lated in Tables 1 through 4. The first test on the unsteered 
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situation is shown in the computer outputs from Figure 5 to 
Figure 24 for three different wind directions and various 
velocities. Tabulated values are shown on Table 1 and 2, 
for yaw and sway. 

The effects are examined by referring to the tables. 

Casel; Head Wind 

As is shown in Tables 1 and 2 and Figures 5-10, effects 
of the head wind were small on the straight forward ship 
motion. In this test, a 5 mile/hr wind velocity is applied 

_3 

to the ship which created -10 rad. variation on the course 

-3 

and drift was -7.4x10 ft. in 5 minuted. 

Case 2: Beam Wind 

Secondly, beam wind disturbances were examined for various 
velocities. According to Tables 1 and 2 and Figures 11 -16, 
the ship greatly deflected from a straight course where the 
wind velocity exceeded t miles/hr. Deviation v/as positive 
22 degree (.395 rad) on the straight course and positive 
9.2 ft drift in 5 minutes when the 7 miles wind velocity was 
applied, as shown in Figures 15 and 16. 

Case 3: Stern Wind 

In this part, the ship was disturbed with stern wind 
which created great difficulty on the straight course. Five 
miles/hr wind velocity created -2.29 degree (-.04 rad) devia- 
tion and t.28 feet drift in 5 minute. But 7 miles/hr wind 
velocity causes the ship to deviate by -332.3 degree from 
the straight course in 4 minute 10 second. Drift was -110 ft 
in 4 minute but this value was reduced to -15 ft in 5 minute. 
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Figure Wind Direction (degree) Wind Velocity (miles) Time (sec) Yaw (rad. 
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Figure Wind Direction (degree) Wind Velocity (miles) Time (sec) Sway (ft) 
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Briefly the system was unstable and the ship was turned 
around. There is no chance to control the ship. Drift was 
negative 110 ft in 4 minute but this value was reduced the 
negative 15 ft in 5 minute. The summarization of this test 
can be stated as the following: computer program (II) re- 

sponses were tabulated for 3 different wind directions and 
various velocities. According to the table, wind distur- 
bance was not so important on the course from any direction 
until 5 miles/hr wind velocity was reached. Critical points 
were 5 and 7 miles/hr wind velocities for stern and beam 
wind respectively. 

Second test: 15 degree rudder displacement was applied 

under the beam and stern wind disturbances. The main pur- 
pose was to show the effect of the rudder displacement on the 
ship course with wind. As is known, positive rudder displace- 
ment causes negative yaw. All responses are shov/n in Figures 
25-30 and tabulated values are shown in Tables 3 and 4. Ac- 
cording to the Tables 3 and 4, when the ship is effected 5 
miles beam wind velocity and positive 15 degree rudder dis- 
placement. So, ship was deflected -327 degree from the 
original path in 5 minute. Drift was 80 ft in 2.5 minute 
and -100 ft in 5 minute. As mentioned earlier, 5 miles/hr 
beam wind with zero rudder displacement has created 8 x 10 ^ 
rad. deviation from the original path, as shown in Table 1. 

But the same condition with 15® rudder displacement has -327® 
deviation from the original path. In this case, the ship 
was so effected with rudder force and started to turn around 
the path. 
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Figure 



Table 3 



Figure 



Table 4 



WD.DI . (degree) 


WD.VEL. (miles) 


Time (sec) 


Yaw (rad) 


180 


3 


300 


-4 


180 


5 


300 


-18 


90 


5 


300 


-5.7 



. Variation of Yaw for Manually Steered Ship (S=15‘’). 



WD.DI . (degree) WD.VEL. (miles) Time (sec) 



Sway (ft) 



180 



3 



85 

240 



{ 



9 

-18 



180 



5 nonlinear oscillation 



90 



r 1 5 0 r 8 0 

^300 -100 



. Variation of Sway for Manually Steered Ship (S=15®) 
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The other situation is a stern wind disturbance with 15® 



rudder displacement. Also known from Table 1, 5 miles/hr 
stern wind velocity with zero rudder displacement has created 
-.04 rad. deviation from the original path. In this case, 

15 degree rudder displacement created -18 rad. deviation in 
5 minute. Ship was turned 360° around the course in 2 minute 
and drifted with divergent oscillation. 

All these tests were made with ship speed being 1 mile/hr 
in the calm water and values of wind velocities were suddenly 
applied. 

B. AUTOMATICALLY STEERED SHIP 

Earlier tests show that the ship was effected by wind 
and rudder. 

In this section, another test was developed to minimize 
the deviation from the original path. For this reason 
rudder deflection is used and is written as a function of t^. 
The equation for the rudder deflection is: 

6 = a 5d + 6d = k^ + kij; 

where 

a = rudder response time constant (time lag) 

k = yaw gain constant 

k^ = yaw rate gain constant 

Thus, the system would be; 
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These control effectors are the function of some parameter 
that could be measured by sensors, usually a motion parameter. 




Let 




Q _ _ WA.S^+ WC.S+WF ^ *^1 (26) 

5d SSA.S“+SSB.S*+SSE.s 2+SBS.S+SES S+a^ 

H = K + K^S 
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So the characteristic equation is 

SSA.S®+(SSB+a^/SSA) . (SSE+a^SSB-a^ . WA. K^) S ^ 

+ (SBS+a^.SSE-aj^.WC.K^-aj^.WA.K)S^+ (SES+a^^.SBS 
-aj^.WF.K^-aj^.WC.K)S+(a^.SES-a^.WF.K) = 0 (27) 

For simplicity, let B^=SSA, B 2 =SSB+aj^ • SSA etc. 

Bi‘ S®+B2S"+ (B3-B4 .K^) S^+ (B5-B6 *Kt-B7K) 

+ (Bg-B^.K^-B^Q.K)S+(Bj^j^-B^2*K) = 0 (28) 

This characteristic equation was calculated using computer 
program 1. Computer program III was used to get transient 
responses in yaw and sway for the automatically steered ship. 
Each curve of yaw or sway was developed with different param- 
eter values (for K and K^) , for different wind speed and 
direction. 

The problem is to minimize the variation on the ship 
course. For this reason, the parameter plane program is ap- 
plied to the system to determine the best value of yaw gain 
constant and yaw rate gain constant. 

The responses of these are shown from Figure 40 to 
Figure 44 for several situations. The values chosen are 
the minimum values shown on the curves according to refer- 

I 

ence [6] . 

The parameter plane outputs are developed by applying 
equation (27) to program IV. Table 5 shows the best values 
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of the yaw gain constant and yaw rate gain constant for the 
automatically steered ship. 

All responses of the automatically steered ship are 
shown in Figure 31 through 33 and tabulated on Tables 6 and 

7. 

According to the Tables (6 and 7) , the variations on 
the course is minimized for the low wind velocities. In 
the head wind disturbance, course variations were reduced 
to 0 rad. in 5 minute. 

One mile/hr wind velocity response reached 0 rad. in 

20 sec., 3 mile/hr gave 0 rad. response in 65 sec. But 5 

-4 

miles/hr wind velocity response was 1x10 rad. in 5 minute 
For the beam v/ind affection, 3 miles/hr wind velocity 
was applied to the ship. Deviation was 0 rad. in 70 sec. 

A stern wind was also tried, 1 mile/hr wind velocity 
caused a deviation which was corrected to 0 rad. in 65 sec. 

As is seen above, the ship trajectory was kept on the 
original path. This desired course or original path was 
assumed to be due North. Also the same test condition was 
applied. 



41 










I 






Figure WD.DI. (degree) WD .VEL. (miles) Time (sec) Yav/ (rad) 



0 1 

180 1 

0 3 

90 3 

0 5 



10 


r-lxlO”"^ 

-0 


300 


50 


.-IxlO""^ 


300 


15 


r-lxlO ^ 
^-0 


300 


20 


r .85x10" 


300 


30 


rlxlO"^ 

\5x10"4 


300 



Table 6. Variation of Yaw for Automatically Steered Ship. 



Figure 



WD.DI. (degree) 


^-JD.VEL. (miles) 


Time (sec) 


Sway (ft) 


0 


1 


300 


lOxlo""* 


180 


1 


300 


13x10"^ 


0 


3 


300 


11.5x10" 


90 . 


3 


300 


-5xl0"^ 



Table 7 . 



Variation of Sway for Automatically Steered Ship. 
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IV. STABILITY ANALYSIS 



Stability may be defined in the following way. Estab- 
lish an equilibrium situation, disturb the equilibrium con- 
dition with arbitrary disturbances. If the ship returns or 
tends to return to the original straight-line condition of 
equilibrium after a disturbances of the smallest amount 
(infinitesimal disturbance) , it is stable. If it departs 
or has the tendency to depart from the original straight- 
line equilibrium condition, it is unstable. 

The unstable ship will go in to a starboard or port 
turn. 
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stability analysis is performed by solving for the eigen 
values of the system matrix ( 4 x 4 ) which is defined on page 
46 . Reference [ 3 ] contains straight forward algebraic solu- 
tion of the linearized equations. According to reference 
[3] a complete algebraic solution of linearized equation (20) 
is obtained with a differential operator. Solution result 
of any independent variable is defined as: 

a,t 0pt o t 

X.... = C..e ^ + C_.e ^ + + C .e 

3 (t) 13 23 m3 



m 

^j(t) = 



C. . e 

i=l 



a, t 
1 



( 29 ) 



Where are arbitrary constants depending upon the initial 

conditions and .... are the roots of the deter- 

minant of the coefficients in the linearized equations of 
motion. 



.. a are different in value, the 
m 



In general o^, 

oit 02 t o t 

terms C^^e , €2^ / ^ can not negate one another. 

The only way that each term can go to zero with increasing 
time is for each of the exponents to be negative. If a is 
a complex number in the form a=a+ib, the following relation- 
ship hold e*^^=e ^=e^^ (cos bt+i sin bt) and the condition 

for stability requires that the real parts of o^, .... 0^ 

be negative as they are complex numbers (the imaginary 0, 
parts of number indicate the angular frequency of oscillation)’. 
The Routh oriterian can also be applied to determine system 
stability. If we look at equation ( 20 ) again, the automat- 
ically steered ship equation would be: 
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ABC.l^-AAD.t|j-BFA.l{J-BGC.V-BGD.U = AAE . S 



-BGE.(|^-BFDi|i +ABD.V-BGF V-BGG U 

• • 

-CAA.ijj-BFGiij- BGA. V+ABE .U-CAB .U 

• • 

QOQ.JJJ+QQO.ij: = QQQ.S + S 

Where 



AAC.S 

CAC.S 



( 30 ) 



QQO = k^, QOQ = k, 



QQQ = tr = a 

Equation (30) will be written in the following way. 



ABC 




-AAD 




-BFA 




0 




-BGC 




0 




-BGD 




0 




-AAE 


0 


,, 


-BGE 


• 


-BFD 




ABD 


• 


-BGF 




0 


• 


-BGG 




0 


• 


-AAC 








ij;+ 




ii^+ 




v+ 




V+ 




U+ 




U+ 




i+ 




0 




-CAA 




-BFG 




0 




-BGA 




ABE 




-CAB 




0 




-CAC 


0 




QQO 




QOQ 




0 




0 




0 




0 




-QQQ 




-1 



g=o 



(31) 

As is seen from equation (31) there are four independent 

variables namely yaw, sway, surge and the rudder angle. The 

form used in equation (29) is applied. 

m . m . m . m . 

iJJ.e V=Z V.e^i^, U=Z U.e ^ , S=l A.e i (32) 



i=l 



i=l 



i=l ^ 



i=l ^ 



Where V^, and are arbitrary constants. 

Substitution of equation (32) into equation (31) ; 



ABCO . ^-AAD . 0 . -BFA 
1 1 




-BGC 




-BGD 




-AAE 


-BGEO . -BFO 
1 

-CAAO . -BFG 
1 




ABDO . -BGF 
1 

-BGA 




-BGG 




-7UVC 




V.+ 

1 


ABEO . -CAB 
1 


U.+ 

1 


-CAC 


QQOa.+QOQ 




0 




0 




-QQQO^-1 



A.= 0 
1 



( 33 ) 
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Solutions of equation (33) are possible for only values of 
0 ^ which satisfy equation (34) . 



(ABCa^^-AADO^-BFA) 

(-BGEa^-BFD) 

(-CAAa^-BFG) 

(QQOo^+QOQ) 



(-BGC) 

(ABDa^-BGF) 

(-BGA) 

0 



(-BGD) 

(-BGG) 

(ABEa^-CAB) 



(-AAE) 

(-AAC) 

(-CAC) 



0 



0 (-QQQo^-1) 



(34) 



Program I gives the solution of this equation as fifth 

order polynomial. For different situations such as wind 

velocity and direction other polynomials were also found. 

As is pointed out in page , it is important for the real 

root to be negative for system stability. 

This test was examined in two situations: 

Situation 1: constant wind directions with varied 
velocities . 

Situation 2; constant wind velocities with varied 
direction. 

The first situation was divided into 3 sections: 

Section a: constant head wind with several velocities. 

Section b: constant beam wind with several velocities. 

Section c: constant stern wind with several velocities. 

The second situation was divided into 2 section: 

Section a: constant 7- miles/hr wind velocity with 
several directions. 

Section b: constant 9 miles/hr wind velocity with 
several directions. 

In order to examine the stability of an automatically steered 
ship in wind, all these situations are solved with computer 
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program I, tabulated on Tables 8-12 and drawn on the 
CALCOMP plotters, from Figures 45 through 49. 

Situation 1: This test was developed, constant wind 

direction disturbeci the ship. Wind velocities changed 
each time . 

a) Head wind disturbance; 

As is shown in Table 8 and Figure 45, the ship 
would be stable in varied head .wind velocities until 13 
miles/hr wind velocity. After this range, the ship could 
be unstable on the horizontal plane. 

b) Beam wind disturbance: 

In this case, the stability range was shorter 
(from 1 to 8 railes/hr) . After the 8 miles/hr wind velocities 
the ship could be strictly unstable, shown on Table 9 and 
Figure 46. 

c) Stern wind disturbance: , 

Stern wind effect on the stability might be the 
shortest range, as shown on Figure 47 and Table 10. This 
range included from 1 to 6 miles/hr. After 6 miles/hr wind 
velocities, the ship could be unstable. 

Situation 2: In this case, constant wind velocity dis- 

turbed the ship. Wind direction was changed for each test. 

a) In 7 miles/hr wind velocity with varied direction. 
From 0 degree through 180 degree wind directions were tried 
in 7 miles/hr wind velocity. According to Table 11 and 
Figure 48, the ship could be unstable when the wind dis- 
turbed the ship from 30 degree, 105 degree and 180 degree 
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wind direction. In the other wind direction and 7 miles/hr 
wind velocity the ship, trajectory could be stable. 

b) In 9 miles/hr wind velocity with varied wind 
direction. As is shown on Figure 49 and Table 12, 9 miles/hr 
wind velocity with each direction were applied to the system. 
The ship could be stable when the wind direction was head, 

60 degree - 75 degree and 120 degree - 135 degree. For the 
rest of the directions (0° - 180®) the ship trajectory was 
unstable . 

All of these results show that some of the wind direc- 
tions and velocities played a very active role. These direc- 
tions are 30 degree, 105 degree and 180 degree. After the 
9 miles/hr wind velocities, many of the wind directions had 
adverse affects on stability. 
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Table 10. Real Roots of Characteristic Equation. 



WD.DI 



Re 0 ^ 



0° 


-7.21 


30° 


+3.46 


45° 


-3.34 


O 

o 


-4.9 


75° 


-4.92 


90° 


- -2.56 


105° 


+ 3 .16 


120° 


-6.22 


135° 


-6.56 


150° 


-.0113 


165° 


-1.6x10" 


180° 


+ .443 



Re Re 



124 


-.014 


0104 


-.138 


082 


-.0114 


088 


-.035 


145 


-.0143 


148 


-7 . 8x10 


0419 


-.419 


0896 


-.0153 


118 


-.0249 


947 


-.198 


956 


-.155 


0286 


-.0286 



Re 0 . 
4 


Re 0 ^ 


-11.11 


-60.25 


-10.215 


-60.44 


-16.42 


-60.38 


-10.63 


-60.3 


-10.6 


-60.26 


-10.4 


-60.4 


-10.2 


-60.56 


-10.8 


-60.29 


-10.92 


-60.28 


-10.331 


-60.364 


-10.3 


-60.44 


-10.286 


-60.62 



Table 11. Real Roots of Characteristic 

Equation, Wind Velocity = 7 miles. 
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WD.DI. 


Re 


Re Oj 


Re O 3 


Re a . 
4 


Re 

D 


0 ® 


-5.08 


-.138 


-.024 


- 10.68 


-60.28 


30® 


+8.53 


-.19 


-.0238 


-10.16 


-60.1 


45® 


+.1054 


-1.25 


-.0275 


-7.85 


-60.8 


O 

O 


-2.25 


-.073 


-.0732 


-10.4 


-60.36 


75® 


-2.14 


-.193 


-.0224 


-10.42 


-60.33 


90® 


+1.94 


-.027 


-.145 


-10.25 


-60.42 


105® 


+8.14 


-.056 


-.056 


-10.16 


-60.62 


120 ® 


-3.07 


-.07 


-.045 


-10.45 


-60.42 


135® 


-4.31 


-.132 


-.038 


-10.55 


-60.34 


150® 


+ 2.8 


-.0368 


-.17 


-10.24 


-60.4 


165® 


+4.1 


-.0128 


-.146 


- 10.2 


-60.47 


180® 


+ 4.67 


-.054 


-.054 


- 10.21 


-60.55 



Table 12. Real Roots of Characteristic 

Equation, Wind Velocity = 9 miles. 
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V. DISCUSSION OF RESULTS AND CONCLUSIONS 



The procedures for system course control and stability 
analysis on the horizontal plane were described. Applying 
the datas which are taken from some references, to the 
three degree freedom equation gave us course control and 
stability analysis computer outputs. As a result of these 
outputs some comments can be made. First of all the un- 
steered case was examined. It is seen that the stern wind 
effects the straight course of the ship the most. When the 
wind velocity reaches 7 miles/hr this effect becomes ex- 
tremely high. In 5 minutes, the deviation of the course 
is about 327 degree. At the same time the effect of the 
head wind on the straight course compared to the stern and 
beam wind seems to be much less; which is seen from the 
computer outputs . 

To see the effect of the rudder deflection on the 
straight course there should be positive 15 degree applied 
to the rudder for beam and stern wind. Under these con- 
ditions negative. yaw and sway were found. When the ship 
is under stern wind, applying positive 15 degree rudder 
deflection will create extremely negative yaw and sway. 

As it is seen in Figures 29 and 30, when the stern wind 
velocity reaches 5 miles in 50 seconds, 15 degree rudder 
deflection will create minus 145 degree deviation from 
the course. 
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mm 






Computer output results were also obtained for an auto- 
matically steered ship. Parameter plane computer programs 
were applied to the closed loop system for choosing the 
best value of yaw rate gain constant and yaw gain constant. 
According to the automatically steered ship computer out- 
puts, in low wind speeds the course variation is minimized. 
Sway can't be controlled very well but the value of the 
drift has decreased compared to the unsteered case. 

In 'a final study, using the eigen value theorem, system 
stability and the critical wind velocity are found and 
tabulated. Computed results show that 30 degree, 105 degree 
and 180 degree wind direction are important on the ship 
course in X miles/hr wind velocities. Also 9 miles/hr wind 
velocity in many direction were so effected. 

All the studies done above are made for a ship which 
has a speed of 1 mile/hr and limited special full loaded 
ship form. 
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APPENDIX 



Hydrodynainic and aerodynamic coefficients where used 
in equation (5) 
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N^, N 2 » are coefficients of aerodynamic forces 

and moments 

Aerodynamic forces and moment 
X_ 
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I _ 
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2" ^ 
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ea p 2 
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= X. ' U ' = X • cos Ipa 
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«al 



= N'u'v'+N'v' = sin 2ijja+Na2 sin ipa 



= y ' V ' = Y ' sin ipa 

-L 3 . 3 
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Which 



X, ' = ' = -0.015 

1 

y ' = Y ' = -0.056 

1 ai 

N, ' = N ' = -0.0046, N^' = -N ' = -0.017 

1 ai 2 a2 
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Figure 5. Variation of Yaw Versus Time for Wind Direction 
000°, Velocity 1 raile/hr. 
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Figure 6. Variation of Sway Versus Time for Wind Direction 
000®, Velocity 1 mile/hr. 
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Figure 7. Variation of Yaw Versus Time for Wind Direction 
000 f Velocity 3 miles/hr. 
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Figure 8. Variation of Sway Versus Time for Wind Direction 
' 000®, Velocity 3 railes/hr. 
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Figure 10. Variation of Sway Versus Time for Wind Direction 
000°/ Velocity 5 miles/hr. 
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Figure 11. 



Versus Time for Wind Direction 
090®, Velocity 3 miles/hr. 
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Figure 12. Variation of Sway Versus Time for Wind Direction 
090®, Velocity 3 miles/hr. 
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Figure 13. ^^J^^tion^of^Yaw Versus Time for Wind Directic 



Velocity 5 miles/hr 
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Figure 14. 



Variation of Sway Versus Time for Wind Direction 
090°, Velocity 5 miles/hr. 
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Figure 15. Variation of Yaw Versus Time for Wind Direction 
090 , Velocity 7 miles/hr. 
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Figure 16. Variation of Sway Versus Time for Wind Direction 
090°/ Velocity 7 miles/hr. 



69 




i 



AW/A 




T/M£- (Sec) 



Figure 17. 



Versus Time for 
180 , Velocity 1 mile/hr. 



Wind Direction 
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Figure 18. Variation of Sway Versus Time for Wind Direction 
180°, Velocity 1 mile/hr. 
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Figure 19. Variation of Yaw Versus Time for Wind Direction 
180 , Velocity 3 miles/hr. 
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Figure 20. Variation of Sway Versus Time for Wind Direction 
180®, Velocity 3 miles/hr. 
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Figure 21. Variation of Yaw Versus Time for Wind Direction 
180®, Velocity 5 miles/hr. 
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Figure 22. Variation of Sway Versus Time for Wind Direction 
180®, Velocity 5 miles/hr. 
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Figure 23. Variation of Yaw Versus Time for Wind Direction 
180®, Velocity 7 miles/hr. 



76 



5UJAY 



/ 

/ 



I 























i 










' 


-jto 




T> 


400 


ffOV 


i — 

o 

*T* 1 — 




\ i 







T/ME' 




Figure 24. Variation of Sway Versus* Time for Wind Direction 
180°, Velocity 7 miles/hr. 
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Figure 25. Variation of Yaw Versus Time for Wind Direction 

090°, Velocity 5 miles/hr. Rudder Deflection 15°. 
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Figure 26 . Variation of Sway Versus Time for Wind Direction 

090 °, Velocity 5 miles/hr, Rudder Deflection 15 °. 
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Figure 27. Variation of Yaw Versus Time for Wind Direction 

180®, Velocity 3 miles/hr, Rudder Deflection 15°. 
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Figure 28. Variation of Sway Versus Time for Wind Direction 
180®, Velocity 3 miles/hr. Rudder Deflection 15® 
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Figure 29. 



Variation of Yaw Versus Time for Wind Direction 
180* , Velocity 5 miles/hr. Rudder Deflection 15°. 
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Figure 30. Variation of Sway Versus Time for Wind Direction 
180°, Velocity 5 railes/hr. Rudder Deflection 15°. 
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Figure 31. Variation of Yaw Versus Time for Wind Direction 
000°/ Velocity 1 mile/hr, Automatically Steered. 
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Figure 32. Variation of Sway Versus Time for Wind Direction 
000°, Velocity 1 mile/hr. Automatically Steered. 
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Figure 33. Variation of Yaw Versus Time for Wind Direction 
000°, Velocity 3 miles/hr, Automatically Steered 
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Figure 34. Variation of Sway Versus Time for Wind Direction 
000°, Velocity 3 miles/hr, Automatically Steered. 
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Figure 35. Variation of Yaw Versus Time for Wind Direction 
000 , Velocity 5 miles/hr. Automatically Steered 



88 




I 




YAW CraJ) 




T/ME 



Figure 36. Variation of Yaw Versus Time for Wind Direction 

090°, Velocity 3 railes/hr. Automatically Steered. 
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Figure 37. Variation of Sway Versus Time for Wind Direction 
090°, Velocity 3 miles/hr, Automatically Steered. 
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Figure 38. Variation of Yaw Versus Time for Wind Direction 
180®, Velocity 1 mile/hr, Automatically Steered. 
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Figure 39. Variation of Sway Versus Time for Wind Direction 
180®, Velocity 1 mile/hr, Automatically Steered. 
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Figure 40. Yaw Gain Constant Versus Yaw Rate Gain Constant 
for Wind Direction 000°, Velocity 1 mile/hr. 
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Figure 41. Yaw Gain Constant Versus Yaw Rate Gain Constant 
for Wind Direction 000°, Velocity 3 railes/hr. 
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Figure 42. Yaw Gain Constant Versus Yaw Rate Gain Constant 
for Wind Direction 000°, Velocity 5 miles/hr. 
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Figure 43. Yaw Gain Constant Versus Yaw Rate Gain Constant 
for Wind Direction 090°, Velocity 3 miles/hr. 
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Figure 44. Yaw Gain Constant Versus Yaw Rate Gain Constant 
for Wind Direction 180°, Velocity 1 mile/hr. 
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Figure 45. Real Versus Wind Velocity for Head Wind. 
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Figure 46. Real Versus Wind Velocity for Beam Wind. 
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Figure 47. Real Versus Wind Velocity for Stern Wind. 
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Figure 48. Real a Versus Wind Directinon for 7 miles/hr Wind 
Velocity. 
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Figure 49. Real Versus Wind Direction for 9 miles/hr Wind 
, Velocity. 
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